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Abstract 

The asymptotic normality of conditional least squares estimators for the offspring variance 
in critical branching processes with non-homogeneous immigration is established, under mo- 
ment assumptions on both reproduction and immigration. The proofs use martingale tech- 
niques and weak convergence results in Skorokhod spaces. 
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1 Introduction 

Kimmel and Axelrod (2002), Pakes (2003), and Haccou et al. (2005) survey applications of 
branching stochastic models in genetics, molecular biology, and microbiology. Yakovlev and 
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Yanev (2006) point out that in vivo cell kinetics requires stochastic modeling of renewing 
cell populations with non-homogeneous immigration. In this line, Hyrien and Yanev (2010) 
model renewing cell populations where the experimentally observable cells are supplemented 
by unobservable cells, for example stem cells. They analyze the population of terminally 
differentiated ohgodendrocytes of the central nervous system and the population of leukaemia 
cells. In both cases the cell population expands through both division of existing (progenitor) 
cells and differentiation of stem cells. This dynamics belongs to branching processes with 
non-homogeneous immigration. The population's viabihty (Jagers and Harding, 2009) is 
preserved by allowing the immigration distribution to vary in time, increasing to infinity on 
average. For semi-stochastic models where immigration depends on the state of the process 
we refer to Cairns (2009). 

In our case of time-dependent immigration (Rahimov, 1995), {-^n,i}n,i>i and {in}n>i are 
two families of independent, non- negative, and integer valued random variables on a prob- 
ability space {il,^,P). We consider discrete-time branching processes with time-dependent 
immigration defined recursively by: 

i=l 

where the {Xn,i}n,i>i have a common distribution for all n and i, and the sequences {Xn,i}n,i>i 
and {Cn}n>i are independent of each other. X„ j is the number of offspring of the ith indi- 
vidual from the [n — l)th generation and ^„ the number of immigrants (or invaders, Haccou 
et al., 2005) joining the population at the time of birth of the nth generation. Z„ rep- 
resents the nth population size and the independence of {Xn,i}n,i>i from {^n}n>i implies 
the independence of reproduction and immigration. In contrast to usual branching process 
models with immigration, we do not assume {^n}n>i to be identically distributed, but that 
the immigration rate varies from generation to generation. 
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Our goal is to estimate the offspring variance 6^ := VarXi i assuming that the immigra- 
tion mean and variance are Icnown. We estimate 6^, based on observing a single trajectory 
{Zi, Z2, . . . , Zn} as n — )■ 00. For reproduction, we assume 

EXi,i = 1 and EX^^^ < 00. (2) 

For immigration, let := E^„, := Var^„, and 7^ := Var(^„ — be finite for every n 
and varying regularly at infinity functions of n defined by: 

an = n"La{n), = n'^Lp{n), and 7^ = n'L^{n), (3) 

where a, /3, 7 are non-negative and \im.x^oo L(^:){cx) / L(.){x) = 1 for any c > 0. The immi- 
gration mean increases to infinity: 

lim an = 00, (4) 

n—>oo 

and the immigration moments, in addition to Eq. (|3]), satisfy: 

(i) lim ^ = or (ii) hm ^ = and lim — ^ = 0. (5) 

n^co nan n^oo na^ n^oo 7* 

Eq. dS]) holds true if /3 < a + 1 or (/3 < 2a; + 1 and a + /S + 1 < 7). Finally, setting 
Vn '■= i^n — OinY — A^, wc assumc that {rjn\n>i satisfies the Lindeberg condition: 

n 

lim^E(^x(l^.|>^))=0, (6) 

fe=l 

where x(^) denotes the indicator of the event A. 

To construct the conditional least squares estimator (CLSE) for the offspring variance. 
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consider Q{k) the cr-algebra generated by {^j}j=i- Denote Efc(-) := E(-|Q'(A;)) and define: 

Mk := Zk — Ek-iZk = Zk — Zk-i — ak- (7) 

Denote 14 := M^-Ek-iM^ The recurrence Eq. ([1]) yields the stochastic regression equation 

Ml = b^Zk-i + (31 + Vk. (8) 

The error terms Vk form a martingale difference sequence {Vk}k>i with respect to the filtra- 
tion {55(A;)}fc>i, that is, Vk is '^{k — l)-measurable and Efc_iV/c = 0. If ak and /3| are known 
and minimizing with respect to 6^ the sum of squares 

n n 
k=l k=l 

the CLSE for 6^ is: 

n n 

Y,{Mi - pi)Zk.i ((^'^ - ^^-1 - "^)' - Pi) ^^-^ 

hi := '-^^n = — -n • (10) 

fc=l fc=l 

For {Z*}„>o, a simple branching process with homogeneous immigration, assuming that 
both offspring and immigration means are known, the CLSE for the offspring variance is 
given by 

n 

^ Y,{Mi)\zu-z:) 

K* ■■= , (11) 

E(^^l - ^nf 

k=l 

where := n^^ J2k=i ^l-i ^1 defined with the first equality of Eq. ([7]), replacing Zk 
by Z^. Yanev (1976/77) and Yanev and Tchoukova-Dantcheva (1986) pioneered the study of 



6^*; Winnicki (1991) proved limit theorems for 6^* assuming finite fourtli moments for both 
reproduction and homogeneous immigration (surveys in Dion, 1993 and Yanev, 2008). Ma 
and Wang (2010) studied the case when these fourth moments may be infinite. 

We extend Winnicki's (1991) findings in the critical case EXi^i = 1 allowing the immigra- 
tion distribution to vary with time such that its mean a„ increases to infinity at a speed that 
could correspond to a near-critical branching population. We prove a limit theorem, which 
completes the results of Rahimov (2008a) in that the limiting distribution of the CLSE for 
the offspring mean depends on the offspring variance 6^. We establish asymptotical normality 
of the CLSE defined by Eq. i^. 

We shall express b'^ — b"^, suitably normalized, as a sum of martingale differences and 
functionals of {Zk}^=i. Then we study the asymptotic behavior of 6^ — 6^ by applying 
either a limit theorem due to Rahimov (2007) giving the limit constant of some functionals 
of {Zk}^^i (Lemma 1) or a general limit theorem, given as Lemma 2, for sequences of 
martingale differences in continuous time (Jacod and Shiryaev, 2003) to each term of the 
sum. 



Define A. 



■n 



n ■ 



.2 . 




I, and On := nAl{nAl + t^)'^. We have 



lim 6. 




■.ee[0,i]. 



(12) 



n 



n^oo nA^ + 



Equality in distribution is denoted by 



d 



Theorem 1 Assume Eq. ([2])- ([6]) hold true. Then 




(13) 
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where A^(0, cr^) is a normal random variable with zero mean and variance 

= {2a + 3)2 f + (1 - \ ■ (M) 

^ ^ \ 4a + 5 ^ ^ 2a + 3 + 77 ^ ' 

For the critical process with homogeneous immigration, Winnicki (1991) established the weak 
limit of 6^* with rate of convergence n^/^. In Theorem 1 the convergence rate is (6'„n)^/2, 
where lim„^oo^n = & [0, 1] and lim„^oo^n^ = oo, provided lim^^oo ''^^«n/7n = 0- 
Corollary 1 Under the assumptions of Theorem 1, 

(i) If limn^^n^al/-/^ = 0, then Eq. (Il3])-(II4]) hold true with ^ = 0. 

(ii) If lim^^oo n^a^/Tn = oo> then Eq. (Il3])-(II4]) hold true with 9=1. 

Example 1 (Poisson immigration) For {^n}n>i independent Poisson variables with mean 
an = 'nfLain) = o{n) — )■ oo as n — oo, Eq. (|4])-(l6]) and the condition (ii) in Corollary 1 
satisfied, if Eq. ([2]) holds true, then Theorem 1 implies Eq. (TT3l)-( |T4|) with 9 = 1. 
Example 2 (Neyman Type A immigration) If Eq. ([2]) holds true and for {^n}n>i 
independent with Neyman Type A distribution given by Ez«" = exp (A„ [e^-^'-^) - l)), 
\z\ < 1, if \n = n^Lxiji) — i- oo, (A > 0) and ipn = n'^L^{n), {ip > 0), then (Johnson et al., 
1993: 371) the rth factorial moment ^n{r) satisfies /i„(r) ~ A^y^J^ for r > 1 and n — > oo. 
Also {r]n}n>i satisfies Eq. (E]). If < A + < 1/2, then Eq. and (I5])(i) hold true and 
hence Eq. ([I3])-([I4]) with 9 = 1 and a = X + Ifl<A + v5< 3/2, then Eq. @ and ([5])(ii) 
are satisfied, which yields Eq. ( JT3l) -( IT4l) with ^ = and a = X + (p. 

2 Preliminaries 

The proof of the theorem uses auxiliary results given in this section. "D " denotes conver- 
gence or equality in the Skorokhod space Z}(M+,R), "P " probability, and "c? " distribution. 
The first lemma summarizes limit results for functionals of Eq. ([1]). Its proof is similar to 
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that of Corollary 2 in Rahimov (2008b) and is omitted here. 

Lemma 1 For the critical process ([1]) with < oo, lim„_^oo = oo and 

P^ina^) = 0, for any continuous function $ on M+ and any sequence {cn}n>Oj 
varying regularly at infinity with exponent p > 0, for t > we have: 

P I ^P$(^"+l)rfy. (15) 



hm —f^Ck^ (^) 







A necessary and sufficient condition for weak convergence in a Skorokhod space of a 
sequence of martingale differences (Jacod and Shiryaev, 2003: Theorem VIII.2.29; Ispany et 
al., 2006) is: 

Lemma 2 (CLT for martingales) For a sequence of martingale differences }fc>i, 
n > 1, with respect to a filtration {55fc}fc>i, such that for all £ > and t > the Lindeberg 
condition 

[nt] 

hm {{U-rx{\U-\ > e) I ^U) = (16) 

k=l 

holds true. Then 

[nt] 

\im Y,U^ = U{t), (17) 

k=l 

where f/(t) is a continuous Gaussian martingale with mean zero and covariance function 
C{t) if and only if for every t > 0, 

[nt] 

\i^Y.^{{Ulf\^l.,)=C{t). (18) 

fc=l 

We use a tilde to indicate that a random variable ^ is centered around its mean, that is 
C = C - EC- Denote Yk.i := X|_., Sj := Yll=i^k,u and recall that 7]ri ■= {in - "n)^- The 



7 



expansion of the error term Vk plays a key role in our analysis: 



Vk = Mi-Ek-iMi (19) 
= {Zk - Zk-i - - h''Zk-i - I3l 

j=2 J i=l 1=1 



Denoting Vn\t) := H^^ Y^t=i ^i*'*' 1 ^ ^ 3, where = nA\ + r^, we decompose the 
normalized sum of the error terms into three parts: 

[nt] 

Vn{t):=—Y.^^ = ^n^^) + V^'^{t) + V^^^{t), t>0. (20) 

k=l 

We shall show that the asymptotic behaviour of l^(t) as n — )■ oo is governed by Vn^\t), 

(2) (3) 

while the contributions of Vn (t) and Vn (t) are negligible. Define 



(21) 



v{t) := w + (1 - ^)r+i^ , t > 0, 

where W{t) is a standard Wiener process and 9 is the limiting constant in Eq. f lT2|) . 
Proposition 1. If Eq. (I2])-(I5]), then for every t > 



hm V«(t)^V(t). (22) 
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Proof. Denote := 2 ^,=2^ ,,S'j_i. The independence of {X„ j}„ j>i and Lemma 1 yield 



[nt] 



EE 

fc=l 



[nt] 



ni] 



~2 



(23) 



fc=i 



k=\ 



[nt] 



nAl 



r ^fc-i(^fc-i ~ 1) 



fc=i 



[nt] 

fe=l 



7.^ 



p 



26^^ 
2a + 3 



Then V^i^^/Hn satisfies Eq. ([HD with C{t) = 2b^et'^''+^{2a + 3)-^ + (1 - e)t^-^\ We shall 
verify Eq. ( |T6i) . Indeed, 



[nt] 



fc=l 



>^ <Tp$^E,_i((T2 + ^)x(|T, + ^fc|>eifO) 

/ y fc=i 

n \ \ / \ 

" k=l ^ 



nix{\vk\>'-^]+nix[\Tk\> 



(24) 



The independence assumption, Chebyshev inequality, and Lemma 1 imply 



h(n) + Un) 



n ^j^^ V \ / \ 



sHn 



[nt] 



< 



" A:=l 



^ (Varfc_iTfc Varr/fc + Varr/^ Var^.iTfc) 



86X(l-^^n)l^^;c.-i(^fe-i-l) 



fc=i 



AH 



-Ik 



4 0. 



(25) 
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For Isin), referring to the Lindeberg condition in Eq. (j6]), we have 



i.|E...(,i,(fei>£f.) 



< 



I -On 



(26) 



[nt] 

Ee.-i {vix{\v,\>'-Y 

k=l 



4 0. 



It remains to show hm„_j,oo hin) = 0. From Burkholder inequahty, for 5 > and p > 1, 



E 



i=2 



2p 



i=2 



(27) 



where Di p > depends on p only. Using Minkovski inequahty, 



E 



E 

i=2 



V2 02 
^k,j'=>j-l 



' < ^ (e (|X,,,f^| Vip^))' = 4E (E| (28) 



i=2 



i=2 



where, by assumption, := E|Xi^ipP < oo for < p < 2. Similarly, we obtain 



E\S,,i\'P < D2,pE 



< /^2,p&2p(j - 1)^ 



(29) 



where 1^2,? > depends on p only. From Eq. (127|) -( 129|) . 



E 



i=2 



2p 



Li=2 



(30) 



where Dp > depends on p only. Applying Eq. ^ with p = (2 + 5)/2 to Ek-i\Tk\'^+^ , for 
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every e > and 6 > 0: 



[nt] 



h{n) < 



k=l 



(31) 



< 



^SH2+S 2^ 
" fc=l 




2+5 



A 0, 



where we have used Lemma 1 and hm„_j,oo An/Hn = 0. Hence, '/Hn satisfies Eq. (HE]) and 
because V^^^ is a martingale difference with respect to {53fc}fc>i, all assumptions of Lemma 2 
are verified. Lemma 2 implies Eq. f l22|) . 

We derive the limit of the normalized sum in Eq. ( I20ll . 

Proposition 2 If Eq. ([2])-([6]), then for every t > 



Proof. From Eq. (!20l) and Proposition 1, in order to prove Eq. ( l32l) . it remains to establish 
that both Vn'^\t) and Vn^\t) are asymptotically negligible. First we show 



lim Vn{t) = V{t). 



(32) 



lim V^^\t) = 0. 



(33) 
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Observing that liirin^oon An /H^ = and applying Lemma 1, we obtain 



[nt] 

EE, 

k=l 



v; 



(3)' 



Hr, 



X 



(3) 



> e 



< 



^ [nt] 



" k=i 



(34) 



C'n ■'^-'1,1 



y4„ n A,, 



4 0, 



where EF^^^ = EXf^ - 4EX3i - 46^ + 36^ + 3 < oo (due to Eq. <^). Then, Eq. (US]) and 
Eq. ([IS]) with C{t) = hold true. Lemma 2 yields 1/i^^(t) 4 and hence Eq. (l33l). 
As for Ki^''(t), applying Lemma 1, we obtain 



(2)' 



fc=l 



X 



V, 



(2) 



> e 



1 / 



An n(5l 



(2) 



(35) 



k=l 



Ar, 



4 0, 



where, according to Eq. ([S]), lirnn^oo Pn^nl An = 0. Therefore, Eq. (ITBl) and ( ITS]) with 
C(t) = are satisfied. From Lemma 2 again, we have lim„^oo Ki^"* (t) = and hence 



lim V^^\t) = 0. 



(36) 



Eq. (132]) follows from Proposition 1, Eq. (133]) . Eq. ( ]36]1 . and Slutsky theorem 
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3 Proof of Theorem 1 

Recalling Eq. ( ITOl) and the definition of V^, we have 

Z^fc=l Z^k=l ^k-l 

First, we examine the asymptotic behavior as — oo of the numerator in Eq. fl37l) . Theo- 
rem 3.1 in Rahimov (2009) and Proposition 2 for t > imply: 



lim = lim ^ = and lim Vnit) = lim V = V{t), (3 

k=l 



nt] 



where the convergence is on D(R_|_,R_(_) and D(]R+,R), respectively. Both limiting processes 
are continuous, referring to Theorem 2.2 in Kurtz and Protter (1991), so that 

lim fz„(t), V;(t), I Zn{u)dVn{u)] = (e+\V{t), I u'^+UViu)], (39) 
V Jq J \ Jo J 

on D(R+,R+ X R X M). Eq. ( I39l) and the continuous mapping theorem (Billingsley, 1968: 
Theorem 5.5) yield: 



^±nz.-. ^ i:^..(^)(K( 

" " A:=l fc=l ^ ^ ^ ^ 



rl„A„ ^-^ ^ — ^ \n J \ \ n I \n 

•1 

Zn{u) dVa{u) 

1 



- K - (40) 



A / u''+^dV{u). 
Jo 

For the denominator in Eq. fl37|) . applying Lemma 1 with t = 1, we have 

n 

" fc=i 
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Finally, Eq. f HOj) and f HTj) and Slutsky theorem yield: 

lim ^ (61 -b^)= lim ^ttCiV^'z^" ^ (2a + 3) , (42) 

which implies the normality of the limit in Eq. ( ITSll . Ito's formula yields: 



1 rl rl 

a+l 



u''+'dV{u) = V{1)- (a + l) / u''V{u)du= (a + l) / (\/(l) - \/(u))u" (iw. (43) 
Jo Jo 

The variable C := /o (V'(l) - V(n))u"ciu is a normally distributed with zero mean and 

EC'= / / s"rE{{V{l) -V{t)){V{l) -V{s)))dsdt. (44) 



For < s < t < 1, 



we calculate 



Emi)-VmW{l)-W{s))) (45) 
= E((\/(l) - \/(t))^ + (\/(l) - VmV{t) - V{s))) 

2a + 3 ^ ^ ^ ^ ^ ^ ' 



EC^ = ( + (1 - 9) ^ + ^ ) , (46) 

^ (a + l)2V4a + 5^ ^2a + 7 + 3/' ^ ^ 



which, taking into account Eq. fj42l) - fj43l) . implies Eq. f|T4l) and completes the proof. 

4 Conclusion 

Branching processes with time-dependent immigration are encountered in a variety of appli- 
cations to population biology. We study conditional least-squares estimators for the offspring 
variance in the critical case, assuming that the immigration mean increases to infinity over 
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time. Theorem 1 establishes the asymptotic normahty of the proposed estimators with 
convergence rate {Onny^"^, where hm^^oo = ^ € [0,1]. A next question concerns the 
conditional consistency of the estimators with different weights. 
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